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Abstract. We study the structure and the internal colour gradients of clus¬ 
ter galaxies from UV to NIR restframe, in the redshift range z = 0.21 — 0.64. 
Structural parameters (half light radius re, mean surface brightness </r>e and 
Sersic index n) are derived for galaxies in the clusters A 209 at z = 0.21 and 
EIS 0048 at z = 0.64. This data set, together with previous data for the cluster 
AC 118 at z = 0.31, constitutes the first large (N ^ 270) sample of cluster galaxies 
whose internal structure in UV, optical (OPT) and NIR (U-, V- and H-band rest- 
frame) can be investigated up to a look-back time of ~ 6 Gyr (flrn = 0.3, Qa = 0.7 
and Ho = 70 Km s“'^Mpc“'^). Galaxies are classified as spheroids or disks accord¬ 
ing to the shape of the light profile, and the evolution of the two populations are in¬ 
vestigated separately. On average, both spheroids and disks are more concentrated 
at longer wavelengths: the galaxy sizes become smaller from UV to NIR, while 
Sersic indices increase. This trend shows an evolution in disks, where the mean 
ratio of optical to NIR Sersic indices decreases from z = 0.31 to z = 0.64. Colour 
gradients are on average negative at all redshifts and are stronger in disks than in 
spheroids. But while for spheroids both grad(UV — OPT) and grad(OPT — NIR) 
are only weakly dependent on redshift, the optical-NIR gradients of disks become 
significantly smaller at z = 0.64. Colour gradients and central colours are com¬ 
pared with models of metallicity, age, and dust extinction gradients. Metallicity 
turns out to be the primary driver of colour gradients in spheroids, the age gradi¬ 
ent being constrained to be smaller than ~ 25%. For disks, two kinds of models fit 
the present data: (i) age gradients (in the range [30, 50] %) with significant dust 
extinction, and (ii) ‘pure’ dust models, in which the gradients of colour excess are 
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a factor of two higher in EIS 0048 than in the other clusters. Since colour gradi¬ 
ents of disks seem not to correlate significantly with inclination, we argue that age 
gradient models could represent a more likely explanation of the present data, in 
agreement with what expected on the basis of hierarchical merging scenarios. 

Key words. Galaxies: clusters: individual: A 209, AC 118, EIS 0048 - Galaxies: 
photometry - Galaxies: evolution - Galaxies: fundamental parameters {colours, 
colour gradients, effective radii, Sersic indices) 

1. Introduction 

Studies of spectral indices and broadband colours have demonstrated that stellar popu¬ 
lations in galaxies are not homogeneous, but are characterized by a spatial variation in 
their physical properties. 

Most of nearby early-type (E) galaxies have negative colour gradients, becoming 
bluer toward the periphery (e.g. Boroson et al. Franx et al. [T?M1 Peletier et 

al. I199()al I1990bl hereafter PDI90 and PVJ90, respectively; Michard I1999II . The pres¬ 
ence of radial changes in colour may be interpreted in terms of both age and metallicity, 
due to the well known degeneracy between the two properties (Worthey et al. I1996II . 
The most effective way to investigate the origins of colour gradients is to study their 
evolution with look-back time, since age and metallicity models give very different pre¬ 
dictions at increasing redshifts (Kodama & Arimoto ll997ll . By adopting this approach, 
Saglia et al. Ii™il and Tamura & Ohta li™il analyzed the colour profile of Es in dis¬ 
tant clusters by using HST optical photometry. They found that galaxies show moderate 
colour gradients also at intermediate redshifts, and claimed, therefore, that metallicity 
is the primary factor driving the observed gradients. The same conclusion was drawn by 
Tamura et al. I|2000|l and by Hinkley & Im 112001 II for early-type galaxies in the field. 
On the other hand, there are also some indications that the population of spheroids at 
intermediate redshifts is not homogeneous with respect to the radial properties of the 
stellar populations. Recently, Bartholomew et al. ii^nnni showed that cluster K-l-A galax¬ 
ies at intermediate redshifts could have somewhat bluer colour gradients with respect 
to normal early-types, while Menanteau et al. llWTll found positive colour gradients for 
spheroids in the Hubble Deep Fields, concluding that a significant fraction of galaxies 
have experienced recent star formation localized in the center. 

Various attempts have been made to interpret the existence of colour gradients of 
early-type galaxies in the framework of different evolutionary models. The presence 
of metallicity gradients is naturally expected within the monolithic collapse scenarios 

Send offprint requests to: F. La Barbera 

* Based on observations collected at European Southern Observatory. Figures G.l, C.2, and 
C.3 are only available in electronic form at http://www.edpsciences.org 
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('Larson flf)74ll . as a consequence of the later onset of the galactic wind in the inner region 
of the galaxy, but it can also be accommodated within the framework of galaxy for¬ 
mation via hierarchical merging (White 119801 Mihos & Hernquist I1994II . Positive colour 
gradients are not predicted by the monolithic collapse while they can be explained as a 
consequence of the centrally peaked star formation produced in merging remnants (Mihos 
& Hernquist 119941 [T99t)|l . Recently, it has been shown that colour gradients of spheroids 
can also be explained within the chemodynamical model for evolution of elliptical galaxies 
(see Friaga & Terlevich 120011 and references therein). 

The origin of colour gradients is crucial to understand the mechanisms underlying the 
formation and evolution of late-type galaxies. Colour gradients in nearby disks have been 
extensively investigated by de Jong ifTMl using optical and near-infrared broadband 
photometry. He found that disk galaxies have more metal rich and older stellar popula¬ 
tions in the center, and that, therefore, both age and metallicity drive the observed radial 
colour prohles. The presence of age gradients is also supported by the increase of the disk 
scale-lengths measured by the Hq, flux with respect to the optical (Ryder & DoDita [r994ll . 
Further investigations have been performed by Gadotti & dos Anjos iiwm . who studied 
a large sample of late-tape spirals in the held by UV - optical photometry. They found 
that most of the galaxies (85%) have negative or null colour gradients, and suggested 
that colour gradients are more sensitive to age than to metallicity of stellar populations. 
Disks in the cluster environment have been studied with NIR (J - K) data by Moriondo 
et al. cnnu, who found very weak or null colour gradients. 

The presence of age gradients is a robust prediction of hierarchical models, in which 
spiral bulges contain older stars than their associated disks, which form by subsequent 
accretion. Evidence in favour of the presence of age gradients for held spirals at interme¬ 
diate redshifts has been found by Abraham et al. CM)- 

A crucial role in the interpretation of colour gradients can be played by dust absorp¬ 
tion. Although the effects of extinction remain still unclear and substantially unresolved 
(e.g. Tamura et al. 120001 de Jong llOOHl Gadotti & dos Anios l^nHTjl . different studies have 
shown that dust absorption could have an important effect in the internal colour distri¬ 
butions of galaxies (e.g. Wise & Silva llOOHI Silva & Elston [TT)94ll . In particular, for what 
concerns disk dominated galaxies, it has been proved that internal extinction strongly 
affects the colour prohle of late-type spirals (e.g. Peletier et al. nnna. For these reasons, 
in the present work we consider simple extinction gradient models for disk galaxies. 

An effective tool to analyze the internal light distribution in galaxies is the comparison 
of structural parameters (like the half light radius and the Sersic index) at different 
wavelengths. This approach allows large samples of galaxies to be studied with a broad 
photometric coverage. By using ground-based data. La Barbera et al. jinna hereafter 
LBM02) obtained the first large sample (N = 94) of both optical and NIR structural 
parameters for cluster galaxies at z 0.3. They suggested that both age and metallicity 
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could drive the observed colour gradients in the cluster populations. In the present work, 
we re-analyze this subject by using a larger sample (N ~ 270) of cluster galaxies at 
different redshifts (z = 0.21,0.31,0.64) with a broader photometric baseline: from UV to 
NIR restframe. We estimate colour gradients for both disks and spheroids, with the aim 
of following the properties of both cluster populations down to z = 0.64. 

The layout of the paper is the following. In Sect. [3 we describe the samples used for 
the present analysis. The derivation of structural parameters is outlined in Sect. 13 while 
Sect. 0] deals with the comparison of UV - NIR parameters at the different redshifts. 
The distributions of colour gradients are analyzed in Sect. [3 and their evolution with 
redshift is compared with predictions of metallicity, age, and extinction models in Sect. El 
Summary and conclusions are drawn in Sect. 0 In the following, we will assume the 
cosmology = 0.3, Ua = 0.7 and Hq = 70 Km s“^Mpc“^. With these parameters, 
the age of the universe is ~ 13.5 Gyr, and the redshifts z = 0.21, z = 0.31 and z = 0.64 
correspond to look-back times of 2.5,3.5 and 6 Gyr, respectively. 


2. The data 

The data relevant for the present study have been collected at the ESO New Technology 
Telescope and at the ESO Very Large Telescope. They include images in B and R 
bands for the cluster A209 at z = 0.21 (RA 2000 = 01:31:57.5, DEG 2000 = —13:34:35), 
in R and K bands for the cluster AG 118 at z = 0.31 (RA 2000 = 00:14:19.5, 
DEG 2000 = —30:23:19), and in V, I (high resolution, HR) and K bands for the cluster 
EIS0048at z = 0.64 (RA 2000 = 00:48:31.6, DEC 2000 = -29:43:18.0). A 209 and AC 118 
are rich, massive clusters of galaxies, with velocity dispersions of 1394 ±94 km s“^ 
(Mercurio et af. 12003^ hereafter MGB03) and 1950 ± 334 km s“^ (Wu et af. n999|| l. 
respectively, and with high X-ray luminosity (Lx 13.8 10^"* erg s“^, see Ebeling et 
al. 119961 and Lx = 62.44 ± 14.41 lO'^"^ ergs s“^, see Wu et al. 119991 respectively). Their 
velocity field and spatial structure suggest the presence of substructures, with a complex 
dynamical scenario (see MGB03 and Andreon l2001l respectively). EIS 0048 is a cluster of 
galaxies having a velocity dispersion of ^ 720 km s“^ (Lobo et al. l2002ll . which is typical 
for clusters at the same redshift, and a spatial structure which shows evident signs of 
substructures (see La Barbera et al. l2003|l . A ROSAT PSPC image, which serendipitously 
contains the cluster field, indicates the presence of X-ray emission. 

The reduction of the photometric data of A 209 and AG 118, and that of the V- and 
K-band images for EIS 0048 is described in Mercurio et al. (2003b, MMM03), Busarello 
et al. (2002) and Andreon (2001), and La Barbera et al. (2003, hereafter LMI03), respec¬ 
tively. The reduction of the I(HR)-band image for EIS 0048 is described in Appendix 
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Table 1. Relevant information on the samples. Column 1: cluster identification. Column 
2: redshift. Column 3: wavebands. Column 4: instruments used for the observation. 
Column 5: pixel scale relative to each waveband in arcsec/pixel. Column 6: FWHM (arc- 
sec) of the seeing for each waveband. Column 7: number of spheroids. Column 8: number 
of disks. For EIS 0048 two values are reported, since the sample with K-band data covers 
a smaller area. Column 9: completeness magnitudes. Column 10: selection used to assign 
cluster membership, Zp and Zg denote photometric and spectroscopic redshifts, respec¬ 
tively. Column 11: references for details on data reduction and on the definition of clus¬ 
ter membership. The following abbreviations are used: Mercuric et al. Il2()()3al MGB03), 
Mercuric et al. 112003 bl MMM03), Busarello et al. BML02), Andreon (2001, AOl), 

La Barbera et al. II2()()3[ LMI03). Note that one out of the N = 20 disks with OPT and 
NIR data in EIS 0048 does not have UV structural parameters, and therefore the total 
number of disks at z = 0.64 is 38. 


CLUSTER 

(1) 

z 

(2) 

Wavebands 

(3) 

Instr. 

(4) 

Scale 

(5) 

Seeing 

(6) 

Nsph 

(7) 

Nd 

(8) 

Me 

(9) 

Sel. 

(10) 

Ref. 

(11) 

A 209 

0.21 

B,R 

EMMI 

0.27,0.27 

0.8, 0.9 

64 

26 

R = 20 

Zs 

MGB03, MMM03 

AC 118 

0.31 

R, K 

EMMI, SOFI 

0.27,0.29 

1.0, 0.8 

76 

15 

R ^ 21 

Zp 

BML02, AOl 

EIS 0048 

0.64 

V,I,K 

FORS2, ISAAC 

0.2,0.1,0.15 

0.6,0.35,0.45 

54,32 

37, 20 

I ^ 22 

Zp 

LMI03 


A. Structural parameters have been derived for galaxies of AC 118 in LBM02^, while 
the surface photometry for A 209 and EIS 0048 is presented in the next section. The 
relevant information on the samples are summarized in Tabled 

In Fig.^ we show the response curves of the optical filters in the restframe for each 
cluster. We see that the R band for A 209 and AC 118, and the I band for EIS 0048 cover 
approximately the same restframe spectral interval, with a central wavelength ranging 

O O 

between ~ 4800 A and ^ 5200 A. These wavebands are close to the V-band restframe 
and will be indicated in the following as optical (OPT) wavebands. The B filter at 

O 

z = 0.21 and the V filter at z = 0.64 have central wavelengths in the range 3200 - 3600A, 
and correspond therefore to U restframe. In the following, they will be indicated as UV 
wavebands. For what concerns the K band (hereafter NIR), it covers approximately H 
and J restframe at z = 0.31 and at z = 0.64, respectively. We note that colour gradients 
and variations in the structural parameters between J and K bands are expected to 
be very small (see Silva & Elston [T994II and therefore the differences of the restframe 
corresponding to the K bands at z = 0.31 and z = 0.64 are negligible for the present 
analysis. The data of A 209 and AC 118 allow UV - OPT and OPT - NIR properties 

^ With respect to LBM02, we excluded two galaxies whose light profiles turned out to be 
strongly affected by nearby objects. 
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A (nm) 


Fig. 1. Restframe filter response of the optical photometric wavebands for A 209, AC 118, 
and EIS 0048. The curves are normalized to 1 at the maximum and are drawn by different 
types of lines as shown in the upper - right. 

of cluster galaxies at 0.2 < z < 0.3 to be studied, while the same wavelength range is 
covered at z = 0.64 by V, I and K bands. Since K corrections do not affect differences 
in colour, the UV-OPT and OPT-NIR colour gradients correspond closely to U-V and 
V-H gradients restframe. We point out that the present samples constitute the first large 
data set of cluster galaxies with structural parameters measured from UV to NIR at 
intermediate redshifts. 

For each cluster we study separately the properties of disks and spheroids. The separa¬ 
tion between the two classes is performed on the basis of the Sersic index n in the optical, 
since the optical parameters are common to each cluster sample. We define as spheroids 
the objects with n > 2, and as disks the remaining galaxies. Following van Dokkum et 
al. this criterion corresponds to distinguishing objects with a low bulge fraction 

(< 20%) from galaxies with a more prominent bulge component. We note that the mean 
value^ of n for disks is 1.10 ± 0.06. By using numerical simulations, we estimated that 
this value describes galaxies with a very low bulge fraction (Lb/Ltot 0.07). Therefore, 
although central bulges can affect the properties of some disks, their effect on the bulk 
properties of disks can be neglected. We also note that, taking into account luminosity 
distances and K-l-E corrections, the completeness magnitudes in Tabled correspond ap¬ 
proximately to the same absolute magnitude limits^: My = —19.70 ± 0.06 at z = 0.21, 

^ The mean value is computed for AC 118 and EIS 0048 galaxies, for which we obtained 
n = 1.2 ±0.1 and n = 1.0 ± 0.05, respectively. For A 209 we found a higher value, n = 1.5 ±0.1, 
due to the presence of a larger fraction of galaxies with n ^ 2. We found, however, that consid¬ 
ering as disks only the galaxies with n < 1.3 in A 209 does not change the estimates (e.g. mean 
colour gradients) given throughout the paper. 

^ Colour gradients of nearby galaxies do not correlate with luminosity (see PDI90), and, 
therefore, slightly different magnitude limits should not affect studies of colour gradients. 








F. La Barbera et al.: Evolution of UV - NIR Structural Properties 


7 


Table 2. Comparison of structural parameters of disks. 



UV- 

OPT 

OPT-NIR 

z 

A(log re) 

A(logn) 

A(log re) 

A(log n) 


mean std 

mean std 

mean std 

mean std 

0.21 

0.01 ± 0.04 0.08 ±0.02 

-0.09 ±0.02 0.10 ±0.04 



0.31 



0.09 ±0.02 0.10 ±0.02 

-0.17 ±0.04 0.20 ±0.05 

0.64 

0.045 ±0.016 0.09 ±0.02 

-0.06 ±0.03 0.15 ±0.02 

0.10 ±0.02 0.13 ±0.07 

-0.34 ±0.05 0.25 ±0.06 


Table 3. Comparison of structural parameters of spheroids. 



UV- 

OPT 

OPT- 

-NIR 

z 

A (log re) 

A(logn) 

A(log re) 

A(logn) 


mean std 

mean std 

mean std 

mean std 

0.21 

0.050 ±0.025 0.12 ±0.02 

-0.02 ±0.07 0.19 ±0.03 



0.31 



0.21 ±0.05 0.26 ±0.04 

-0.03 ±0.02 0.14 ±0.04 

0.64 

0.12 ± 0.03 0.24 ±0.03 

0.05 ±0.06 0.28 ±0.04 

0.080 ±0.015 0.14 ±0.04 

-0.19 ±0.05 0.20 ±0.20 


My = —19.7 ± 0.1 at z = 0.31, and My ^ —20.0 ± 0.2 at z = 0.64. K and evolution cor¬ 
rections were obtained by the GISSELOO synthesis code (Bruzual & Chariot fl993ll . with 
spectral models having an exponential star formation rate, and the uncertainties on My 
were estimated by computing E corrections for different models, exploring a wide range 
of stellar population parameters (i.e. metallicity, formation epoch, time scale of star for¬ 
mation). 

3. Derivation of structural parameters 

Structural parameters were derived for A 209 and EIS 0048 by following the procedure 
described in LBM02. We give here only a brief outline of the method, while details can 
be found in that paper. 

Galaxy images are fitted by the convolution of an accurate model of the Point Spread 
Function (PSF) with a model of the galaxy light distribution. In this work, we consider 
Sersic models: 


l(r) = lo • exp(-b • (r/re)^/"). 


( 1 ) 
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where r is the equivalent radius, re is the effective radius, Iq is the central surface bright¬ 
ness, n is the Sersic index, and b is a constant (b ~ 2n — 1/3, see Caon et al. ll9H3ll . The 
fits were performed by constructing an interactive mask for each galaxy, and by mod¬ 
elling simultaneously overlapping objects. The PSFs were modelled by a multi-Gaussian 
expansion, taking into account possible variations with the chip position. With respect to 
LBM02, we had also to consider that the PSF of the FORS2 I-band images showed sig¬ 
nificant deviations from the circular shape. Details can be found in Appendix B, where 
we illustrate the PSF modelling for the I-band images of EIS 0048. The htted surface 
brightness profiles are shown in Appendix C. 

For AC 118, the typical uncertainties on the structural parameters were estimated 
in LBM02 by comparing the ground-based parameters with measurements obtained by 
HST data. The quoted uncertainties amount to ~ 38% in I'e and 20% in n. For 
A 209, the errors on re and n were estimated by using numerical simulations. We found 
(5(re) ~ 20% and 5(n) ~ 13%. For EIS 0048, due to the very good quality of the images, 
the uncertainties on re and n turned out to be smaller. By comparing measurements 
obtained for galaxies in common between different pointings, we obtained br^lve ~ 15% 
and 6n/n ~ 10%. The estimates of (5(re) and <5(n) were used to derive the uncertainties 
on the colour gradients (see Sect. 0. 

4. Structure of galaxies from UV to NIR 

The shape of the light profile for the Sersic models can be fully characterized in terms of 
the effective radius and the Sersic index. The differences of re and n between the various 
wavebands are shown in Figs. □ [HQ and|3 Differences are always computed between 
shorter and longer wavelengths. Figs. [21 and 0] show the UV - OPT comparison for 
A 209 and EIS 0048, respectively, while Figs. 0 and [31 show the OPT - NIR parameters 
for AC 118 and EIS 0048, respectively. In each figure, left and right panels refer to disks 
and spheroids, while upper and lower panels refer to re and n, respectively. In order to 
analyze the overall properties of each distribution, we derived the corresponding mean 
value and standard deviation by applying the bi-weight statistics (e.g. Beers et al. ITMTIi . 
that has the advantage to minimize the effect of outliers in the distributions. The values 
and the corresponding uncertainties, obtained by the bootstrap method, are shown in 
Tables 121 and 121 for disks and spheroids, respectively. 

4.1. Disks 

By looking at Table (21 we see that disk galaxies tend to have more flat and extended 
profiles in the UV than in the optical. We find, in fact, n'^'^ < n®^"^ and r^'^ > r®^"^ 
both for A 209 and for EIS 0048, although the mean value of A(log re) at z = 0.21 is not 
significantly different from zero within the corresponding uncertainty. The mean value of 
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Log(r.“7r,‘>'=b Log(r.-Vr.°'=b 




Fig. 2. Comparison of UV and optical structural parameters of spheroids (left panels) and 
disks (right panels) for the cluster A 209. Upper and lower panels refer to the effective radii 
and Sersic indices, respectively. The vertical grey line marks the mean of the Gaussian 
fit (dot-dashed line) to the distribution. 
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Fig. 3. The same of Fig. for optical and NIR structural parameters of galaxies in 
AC 118. 
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Fig. 4. Same of Fig. [21 for the cluster EIS 0048. 
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Fig. 5. Same of Fig. Ofor the cluster EIS 0048. 

nUV/nOPT jg ^ Q_g 2 ^ Q_ 2 i and ~ 0.87 at z = 0.64, while the value of r^v/rOPT ranges 
from ^ 1.02 for A 209 to ~ 1.11 for EIS 0048. We note that these difference are within 
Icr and, therefore, are not very significant. At z = 0.64, disks have larger dispersion in the 
log(n^^/n®^'^) diagram, with a standard deviation that is ^ 3/2 of that at z = 0.21. The 
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effect is clearly seen by comparing Figs. Eland01 and is particularly significant in view of 
the fact that the uncertainties on the structural parameters of EIS 0048 are smaller with 
respect to those of AC 118 and A 209. 

Disk galaxies in AC 118 and EIS 0048 have, on average, r®^"^ > r^^^ and 
nOPT < nNiR (see Table EJ, and show, therefore, the same trend observed in the UV 
- OPT comparison: at larger wavelengths, the light profile becomes more concentrated. 
On average, the value of r®P'P/r^^P is ~ 1.23 at z = 0.21 and ~ 1.26 at z = 0.64, and 
nOPT/nNiR decreases from ~ 0.7 to ^ 0.45. While the behaviour of i-Opt^j-NIR similar 
to that found in the UV - OPT, we see that the value of n°P'P/n'^^P shows a much larger 
variation, with a strong decrease for EIS 0048. We will further investigate this point in 
the analysis of colour gradients fSect. IO~Tll . It is also interesting to note that the distri¬ 
bution of Alog(n) at z = 0.31 has a significant tail toward negative values, showing that 
also for AC 118 there are some disks whose NIR light is very peaked in the center with 
respect to the optical. 

Finally, we point out that, since galaxies of AC 118 and EIS 0048 are selected by 
the photometric redshift technique, we cannot exclude that the previous results can be 
affected by the contamination of the cluster samples at z = 0.64 by field galaxies, whose 
properties, also at low redshift, are known to be more disperse and heterogeneous with 
respect to those of galaxies in the cluster environment. The role of field contamination 
will be further addressed in Sect. IHTuI in the comparison of colour gradients. 

4.2. Spheroids 

By comparing UV and optical parameters we see that, on average, spheroids have 
r^^ > The mean value of r^^/r®^"'" is ~ 1.12 at z = 0.21 and ~ 1.32 at z = 0.64, 

while nP^/n^P"'" is consistent with zero at both redshifts within the corresponding un¬ 
certainty. These results indicate that the value of r^^/r^P'^ could increase slightly with z 
while possible variations of nP^/n^P"”" are too small to be significant. It is also interesting 
to note that the dispersions of the two distributions are remarkably larger at z = 0.64, 
ranging from 0.12 to ~ 0.24 for log re and from ~ 0.19 to ~ 0.28 for logn. This is 
similar to what found for the disks and implies that the structural properties of cluster 
galaxies become more heterogeneous at z = 0.64. Since a low fraction of spheroids at 
z = 0.64 are expected to be field galaxies 10 ~ 15%, see Sect. 8.1 of LMI03), these 
results should not to be affected by field contamination. 

Eor what concerns the OPT - NIR parameters, we find that, on average, r^P"”" > r^^^ 
and nOP^ < n^^^ for both AC 118 and EIS 0048. The mean value of nOP^/n^™ varies 
from ^ 0.93 at z = 0.31 to ^ 0.65 at z = 0.64, while r®P'P/r^^P seems to be larger at 
z = 0.31, ranging from ~ 1.20 for EIS 0048 to ~ 1.6 for AC 118. This trend is opposite 
to that found for disks and for spheroids in the UV - OPT. As it is evident from the 
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Table 4. UV — OPT and OPT — NIR colour gradients and central colours of disks. 


grad(UV-OPT) | grad(OPT-NIR) | UV - OPT | OPT - NIR 
mean std mean std mean mean 


0.21 -0.4 ±0.1 0.40 ±0.05 . . 2.0 ± 0.3 

0.31 . . -0.68 ±0.12 0.87 ±0.30 . 3.80 ±0.13 

0.64 -0.46 ±0.12 0.55 ±0.06 -1.2 ±0.1 0.5 ± 0.3 2.2 ± 0.4 4.0 ± 0.4 


Table 5. UV — OPT and OPT — NIR colour gradients and central colours of spheroids. 


grad(UV-OPT) | grad(OPT-NIR) | UV - OPT | OPT - NIR 
mean std mean std mean mean 


0.21 -0.17 ±0.04 0.35 ±0.07 . . 2.4 ±0.1 

0.31 . . -0.40 ±0.14 0.6 ±0.1 . 3.8 ± 0.2 

0.64 -0.33 ±0.17 0.80 ±0.14 -0.38 ±0.15 0.56 ±0.15 2.6 ± 0.2 3.6 ± 0.2 


large dispersion of the distribution in Fig.|31(upper - right panel), this effect could be not 
real but due to the uncertainties on the estimated structural parameters at z ~ 0.3. We 
note, in fact, that also the standard deviation of the log(r®^"'"/r^^^) distribution shows 
a different trend with z with respect to that found for the disks and for the spheroids in 
the UV - OPT: it seems to decrease from z 0.31 to z ~ 0.64. More NIR data at z < 0.6 
would be needed in order to further investigate the previous trends. 

5. Colour gradients 

The dependence of structural parameters on the waveband carries information on the 
properties of stellar populations at different radii from the center of the galaxies. To 
investigate this subject, we derived the internal UV - OPT and OPT - NIR colour 
gradients for A 209, AC 118 and EIS 0048. We used the procedure described in LBM02: 
structural parameters were used to derive the colour fi 2 — Mi of the galaxy at an inner 
radius r, and at an outer radius ro. Here, the symbols 1 and 2 denote the two wavebands. 
The colour gradient was obtained by the logarithmic gradient of M 2 — Mi (see Eq. (4) in 
LBM02). In order to make a direct comparison with previous works (PDI90, PVJ90), we 
adopted ro = re and ri = 0.1 • rg. 


5.1. Distributions of colour gradients 

Colour gradients for A 209, AC 118 and EIS 0048 are shown in Fig. El Hand|H1 where 
upper and lower panels denote disks and spheroids. The bi-weight statistics was applied 
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Fig. 6. Distribution of UV - OPT colour gradients for A 209. Upper and lower panels 
refer to disks and spheroids, respectively. 


in order to characterize each distribution by the relative mean and standard deviation, 
whose values are summarized in Tables 0] and El for disks and spheroids, respectively. In 
the same tables, we also report the mean values of the UV-OPT and OPT-NIR central 
colours (/i 2 — /ii at r = ly). The relative distributions and the standard deviations are 
not shown for brevity, since they are not relevant for the following analysis. The UV- 
OPT and OPT-NIR central colours will be used in Sect. El to discuss the origin of colour 
gradients. We note that the UV — OPT colour indicates B — R and V — I observed colours 
for A 209 and EIS 0048, respectively, while the OPT — NIR colour denotes R — K and 
I —K for AC 118 and EIS 0048, respectively. As mentioned in Sect. 2, the values of 
grad(UV — OPT) and grad(OPT — NIR) correspond approximately to grad(U — V) and 
grad(V — H) restframe, respectively. 

5.1.1. Disks. 

The mean value of the UV —OPT colour gradients (hereafter grad(UV — OPT)) amounts 
to ~ —0.4 and is fully consistent between A 209 and EIS 0048. This is due to the fact that 
the difference between OPT and UV structural parameters does not change significantly 
with redshift (Sect.^Tj. We note, however, that the distribution of EIS 0048 is more 
disperse and asymmetric with respect to that of A 209, with a significant fraction of 
galaxies 40%) having grad(UV — OPT) > 0. This is a consequence of what found in 
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grad(OPT-NIR) 



Fig. 7. The same of Fig. Elfor the OPT - NIR colour gradients of galaxies in AC 118. 




grad(UV-OPT) grod(OPT-NIR) 




grad(UV-OPT) grad(OPT-NIR) 


Fig. 8. The same of Fig. Elfor the UV-OPT and OPT-NIR colour gradients of galaxies 
in EIS 0048. 

Sect. 14. ll some disks at z = 0.64 have a light prohle more peaked toward the center in the 
UV with respect to the optical, and show, therefore, positive values of grad(UV — OPT). 

The most striking difference in OPT - NIR colour gradients (grad(OPT — NIR)) 
between AC 118 and EIS 0048 is the mean value of the two distributions. As shown in 
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Table 01 the mean value of grad(OPT — NIR) is ~ —0.68 at z = 0.31 and ~ —1.2 at 
z = 0.64, implying that there is a sharp evolution in the OPT - NIR colour gradients 
from z = 0.31 to z = 0.64: the NIR light of disks become much more peaked in the center 
at higher redshift (cfr. Sect. 14.111 . 

In order to account for field contamination^ in the previous results, we proceeded as 
follows. As discussed in Sect. 2 of Merluzzi et al. II2003II . the effect of contaminants for 
AC 118 is expected to be negligible, while for EIS 0048 the fraction of field contaminants 
can be estimated as described in LMI03. By considering the colour distribution of disks 
at z = 0.64, we estimated that ^ 6 (out of 20) galaxies in the OPT+NIR sample are 
expected to be field galaxies, and applied a statistically subtraction procedure to correct 
the distributions of colour gradients. The mean value of grad(OPT — NIR) corrected for 
field contamination amounts to —1.35 ± 0.07, which is in agreement with that given in 
Table 0 ] The same result holds for UV-OPT colour gradients. 

5.1.2. Spheroids. 

By looking at Table El we see that the mean value of grad(UV — OPT) is ~ —0.17 
at z = 0.21 and ^ —0.33 at z = 64. While this difference is not particularly significant 
(< Ict), the difference in the scatter of the two distributions is much more evident, as 
it is clearly seen by comparing Fig.E](lower panel) and Fig. 0] (lower - left panel). The 
standard deviation of grad(UV — OPT) for FIS 0048 is more than double with respect 
to A 209. This is due to the fact that, as discussed in Sect. 021 the dispersions of the 
A(logre) and A(logn) distributions increase significantly from z = 0.21 to z = 0.64. 

The mean values of grad(OPT — NIR) are ~ —0.51 for AC 118 and ~ —0.38 for 
FIS 0048, and are consistent within the relative uncertainties. We conclude, therefore, 
that no strong evolution in the OPT - NIR colour gradients is found between z = 0.31 
and z = 0.64 for the population of spheroids. 

6. Evolution of colour gradients 

To investigate the origins of the observed gradients, we compare our results (colour 
gradients and central colours) with the predictions of population synthesis models. 

6.1. Models. 

Colour gradients are modelled by using the GISSFLOO synthesis code (Bruzual & Chariot 
Galaxies are described by an inner stellar population, with age® Ti and metal- 
We point out that, since galaxies were selected by the photometric redshifts, field contam¬ 
inants in the present sample are expected to lie approximately at the same redshift of cluster 
galaxies. 

® Here and in the following age refers to z=0. 
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licity Zi, and by an outer stellar population, with age To and metallicity Zq. Each of 
them is characterized by an exponential star formation rate e”*/’’. For spheroids, we use 
r = 1 Gyr, which gives a suitable description of the evolution of their integrated colours 
(e.g. Merluzzi et al. I2()():-{|I . while for the disks we construct models with both r = 1 Gyr 
and r = 3 Gyr, in order to account for a star formation rate that lasts longer in time. 
For disks, we also consider the effect of dust absorption, by describing the inner and 
outer populations with colour excesses E(B — V)i and E(B — V)o, respectively. The ab¬ 
sorption in UV-OPT and OPT-NIR colours is computed by using the extinction law from 
Seaton dza. For A 209, AC 118 and FIS 0048, the colours of the inner and outer pop¬ 
ulations are obtained by using the instrumental response curves®, while colour gradients 
at z ~ 0 are computed by using the U-, V-, and K-band filter curves from Buser dTHli, 
Azusienis & Straizys dSH), and Bessell & Brett dHHI, respectively. 

Each gradient model is characterized by a set of parameters Pk (Zi,Ti,...), whose 
value is derived from the observations. To this aim, we minimize the chi-square: 

x"(Pk)=E[(Ek-Ok)>^], (2) 

where Ok and CTk are the observed quantities (central colours and colour gradients) and 
the relative uncertainties at the different redshifts, while Ek are the model predictions 
for Ok- 


6.1.1. Spheroids. 

We consider three kinds of gradient models: metallicity (Z) models, age (T) models, and 
metallicity-|-age (TZ) models. The relevant information on the models are summarized 
in Table El For models (Z), we assume that the inner and outer populations are coeval, 
with Ti = To = 13 Gyr, while the values of Zi and Zq are changed. We consider two 
cases: model (ZI), whose inner and outer metallicities are constrained by the observa¬ 
tions, and model (Z2), which describes a strong metallicity gradient, with Z, = 3Zq and 
Zo = l/20Zi. We note that inverse metallicity models (Zq > Z,) are not considered here 
since they are not consistent with the data (see Sect. 1^ . For models (T), we assume 
that the inner and the outer stellar populations have the same metallicity, Zi = Zo, while 
To and Ti are changed. We consider two cases: model (Tl) has Ti = 13 Gyr, Zi = Zq and 
To < Ti, while model (T2) has Zi = 3Zq, To = 13 Gyr and a younger population in the 
inner region. The value of To for model (Tl) is derived on the basis of the observations 
while for model (T2) we chose Ti = 0.8 ■ To. We note that model (T2) describes a galaxy 
in which the stars redden toward the outskirts, and is considered in order to investigate 
the origin of positive colour gradients. Finally, in model (TZ) the inner population is old 
(Ti = 13 Gyr), while Zi, Zo and To are obtained from the observations. 

® The response curves were retrieved from the ESO Exposure Time Calculators at 
http: / / www.eso.org/observing/etc / 
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Table 6. Parameters of colour gradient models of spheroids. Column 1: model. Columns 
2, 3: inner and outer metallicities. Columns 4, 5: inner and outer ages. Notes: (a) con¬ 
strained by both central colours and colour gradients, (b) the parameters are chosen in 
order to describe extremely negative colour gradients, (c) constrained by grad(UV-OPT) 
at z = 0, (d) the parameters are chosen in order to describe positive colour gradients. 


MOD 

Zi 

Zo 

Ti 

To 

(1) 

(2) 

(3) 

(4) 

(5) 

Zl 

3/2Z0(“^ 

o.eZiU) 

13 Gyr 

13 Gyr 


32© 

l/20Zi 

13 Gyr 

13 Gyr 

Tl 

Z© 

z© 

13 Gyr 

0.62TiF) 

'P2U) 

3Zq 

32© 

O.STi 

13 Gyr 

TZ 

3/2Z©(“^ 

1/2Z©(“^ 

13 Gyr 

0.85TiU) 


6.1.2. Disks. 

The relevant information on the gradient models are summarized in Table[7| We consider 
the following models: metallicity (Z) and age (T) models, a ‘protracted’ age (Tl) model, 
a protracted age-bdust (TD) model, and a ‘changing’ dust (D) model. For (Z) and (T) 
models, we consider only the models (Zl) and (Tl) introduced for the spheroids, since 
(Z2) and (T2) are not relevant for the present discussion. Model (Tl) is similar to (Tl) but 
is characterized by a more protracted star formation (r = 3 Gyr) in both the inner and 
outer regions. The same value of r is used for models (TD) and (D). Model (TD) differs 
from (Tl) for the presence of dust in the inner and outer regions. The values of E(B — V), 
and E(B — V)o are assumed to be equal and are derived from the observations. Einally, 
model (D) has Z, = Zq = Zq, T, = Tq = 13 Gyr, and E(B — V)o = 0. The variation of 
colour gradients with redshift is described by changing the amount of dust absorption in 
the inner region between z < 0.31 and z = 0.64. 

6.2. Comparison with observations. 

The observed colour gradients are compared with the models in Figs. El and rm for 
spheroids and disks, respectively. The curves in the figures were obtained by interpolating 
the values of the colour gradient models at z = 0, 0.21, 0.31 and 0.64. 

To further constrain the models, we also considered measurements of colour gradients 
for nearby galaxies. For the spheroids, we use the UV - OPT values for a sample of 39 
nearby ellipticals by PDI90, and the OPT - NIR colour gradients obtained by PVJ90 for 
19 elliptical galaxies. For the disks, an homogeneous sample for which colour gradient 
measurements are available at z ~ 0 both in UV - OPT and in OPT - NIR is not 
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Table 7. Parameters of colour gradient models of disks. Column 1: model. Columns 2, 
3: inner and outer metallicities. Columns 4, 5: inner and outer ages. Columns 6 , 7: inner 
and outer colour excesses. Column 8 : time scale of the star formation rate. Notes: (a) 
constrained by both central colours and grad(UV — OPT) at z = 0.21 and z = 0.64, (b) 
constrained by the colour gradients, (c) constrained by both central colours and UV- 
OPT, OPT-NIR colour gradients, (d) the first value of E(B — V)i refers to z = 0.21 and 
z = 0.31, while the second to z = 0.64. 


MOD 

Zi 

Zo 

Ti 

T 

-L 0 

E(B - V)i 

E(B - V)o 

T 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Z1 

3Z0(“) 

O.8Z0(O 

13 Gyr 

13 Gyr 

0. 

0. 

1 Gyr 

T1 

O.44Z0('’) 

O.44Z0f*') 

13 Gyr 

0.55Tif*’) 

0. 

0. 

1 Gyr 

T1 

O.44Z0(‘') 

O.44Z0f*'^ 

13 Gyr 

0.55Ti(*’) 

0. 

0. 

3 Gyr 

TD 

O.44Z0C) 

O.44Z0C) 

13 Gyr 

0.55TiC) 

O.lsC) 

0.180 

3 Gyr 

dO) 

Z0 

Z0 

13 Gyr 

13 Gyr 

0.15,0.30(‘*> 

0. 

3 Gyr 


available, and therefore we consider, as reference, the values relative to the spheroids. In 
the figures, the black symbols denote the mean values of the colour gradients and the 
relative uncertainties (Icr standard intervals), while the grey bars, which mark the 32% 
percentile intervals of the gradient distributions, describe the range of values of colour 
gradients at each redshift. 

The calibration of the models was performed by minimizing the ‘distance’ (see Eq. ©) 
of the model predictions from the observed quantities (central colours and/or colour 
gradients), as summarized in Tables El and [3 

6.2.1. Spheroids. 

By looking at Fig. El we begin to note that the size of the grey bars at z ~ 0 are 
much smaller than those at higher redshift. Considering the typical uncertainty on the 
estimated colour gradients, we see that the difference is not due to the uncertainty on 
the estimate of structural parameters. The simplest explanation resides in the different 
selection criteria of our samples and those of PDI90 and PVJ90: the samples at z ^ 0 
consist, in fact, of morphologically selected elliptical galaxies, that, in many regards, 
are known to have very homogeneous properties (see e.g. van Dokkum et al. I1998|l . To 
discuss the dispersion in the values of colour gradients, we will consider, therefore, only 
the estimates at z > 0 . 

Fig. El shows that model Z1 gives a good description of all the observed colour gradi¬ 
ents. Moreover this model produces the following central colours: (UV — OPT)g 21 ~ 2.65, 
(UV - OPT)o 04 = 2.58, (OPT - NIR)g 34 ~ 3.7, (OPT - NIR)g 54 = 3.2. These values 






F. La Barbera et al.: Evolution of UV - NIR Structural Properties 


19 
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z 

Fig. 9. UV — OPT and OPT —NIR colour gradients of spheroids as a function of redshift. 
The black and the grey error bars denote, respectively, the uncertainties on the mean 
value of the colour gradient and the 32% percentile interval of the distribution of colour 
gradients. Metallicity and age models are represented by the dashed and the continuous 
curves. The typical uncertainties on individual colour gradients are shown by the error 
bars in the lower - left of the bottom panel for AC 118 and EIS0048 (left and right). 
Typical errors for A 209 are intermediate. 

are consistent (within at most 2a) with the observations (see Table ISJ . We also note 
that model Z1 produces very small OPT - OPT colour gradients in the considered red- 
shift range (e.g. grad(R — I)o. 3 i = —0.04), and is therefore in agreement with optical - 
optical estimates obtained by other studies for cluster galaxies at intermediate redshift 
(e.g. Saglia et al. 120001 Tamura & Ohta l20077ll . The age gradient model (Tl) is consistent 
with the values of grad(OPT — NIR) at z ~ 0 and at z ~ 0.31, as discussed in LBM02. 
However the other values of colour gradients exclude definitively this model, confirming 
that metallicity is the primary cause driving the colour gradients. For what concerns 
model TZ, we see that the combined effect of age and metallicity can explain to some 
extent both the colour gradients and the central colours at all redshifts. This model pro¬ 
duces, in fact, the same central colours of model Z1 and is therefore in agreement with 
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the values in Table |S1 On the other hand, the very rapid decline of grad(UV — OPT) 
and grad(OPT — NIR) produced by age gradients (model Tl) as a function of redshift, 
strongly constrain the maximum allowed age variation, A(T)niax = (T, — To)max, from 
the galaxy center to the outskirts. In order to estimate A(T)max, we shifted the values of 
the colour gradients and of the central colours according to the relative uncertainties by 
using normal random deviates. For each iteration, we re-derived the values of Zo and To 
of model (TZ). We found that the values of Zo and To change according to the relation 
(5(logTo) ~ —2/3<5(logZo), which describes the well known age-metallicity degeneracy 
(Worthey et al. I199()ll . and that ~ 90% of the To values are greater than 10 Gyr, 
constraining A(T) to be smaller than ^ 23% of Tj. 

For what concerns the grey bars, we see that strong metallicity gradients (model Z2) 
are able to describe the tail toward negative values of the colour gradient distributions. 
The OPT - NIR diagram excludes, however, metallicity gradients stronger than that 
of model Z2, setting the constraint Zo>0.15-Zq. Positive values of grad(UV —OPT) 
can be described by age models with T, < To, like model T2. Inverse metallicity 
models (Z, < Zq) give positive values of grad(UV — OPT) but too high values of 
grad(OPT — NIR), and cannot fit, therefore, both gradients simultaneously. The OPT 
- NIR diagram excludes colour gradients significantly higher than those of model T2, 
setting the constraint Ti > 0.8 • To. 

6.2.2. Disks. 

By looking at Fig.^1 it is evident that metallicity models fail to reproduce the variation 
of grad(OPT — NIR) with redshift. Moreover, metallicity gradients in local disks are 
known to be small (e.g. Zaritsky et al. cnniii and therefore they cannot account for the 
large colour gradients observed in the present samples. 

On the contrary, age models (Tl and Tl) give a better description of the data. In 
model Tl, grad(UV — OPT) varies too rapidly with redshift, deviating by ^ 3a from the 
observed gradient a,t z = 0.64. Increasing the r value, the variation of grad(UV — OPT) 
becomes milder, and the agreement improves^. Model (Tl) is able, in fact, to repro¬ 
duce within at most 1.3a the values of grad(UV — OPT) and grad(OPT — NIR) at any 
redshift, and gives colour gradients at z ~ 0 which are only slightly bluer than those of 
the spheroids. Although age models reproduce the observed colour gradients, they fail 
to describe the red central colours of disks (Table 0J, which are not significantly dif¬ 
ferent from those of the spheroids (Table 0. Model Tl gives (UV — OPT)p 21 ~ 2.0, 

^ The UV-OPT colour is more sensitive with respect to OPT-NIR colours to the presence of 
young stars. For high values of r, the inner and outer stellar populations protract the formation 
of fresh stars with redshift, and therefore the evolution of grad(UV — OPT) becomes smaller 
with respect to that of grad(OPT — NIR). 
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Fig. 10. Same of Fig. 13 for the disks. 

(UV - OPT)o 54 ~ 1.9, (OPT - NIR )5 34 ~ 3.0, (OPT - NIR )5 54 = 2.5, which are not 
consistent with the observations. We also note that, by increasing the inner and outer 
metallicities, the central colours are better reproduced, but the description of colour 
gradients worsens®. By fitting Ti, To, Z,, Zq to both the colour gradients and the central 
colours, we found that these parameters are not able by themselves to describe all the 
observed quantities at better than 2.5a. We conclude, therefore, that stellar population 
parameters cannot explain alone our data. 

The possible solution resides in dust absorption, which is known to produce 
strong colour gradients in nearby late-type spirals (Peletier et al. I1995II . Model TD 
gives the same observed gradients of model Tl, and produces the following central 
colours: (UV - OPT )5 21 ~ 2.4, (UV - OPT )5 54 ~ 2.3, (OPT - NIR )5 34 ~ 3.9, 
(OPT —NIR)Qg 4 = 3.3, in good agreement (< 1.5 (t) with the values in Tabled We 
point out that the present data can also be explained by models in which both age 
and dust content vary toward the periphery. In order to investigate this point, we 
® In fact, grad(UV — OPT) is more sensitive to the metallicity with respect to 
grad(OPT — NIR), in the sense that increasing Zi and Zo, the UV-OPT gradients decrease 
more rapidly than grad(OPT — NIR). 
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considered To and E(B — V)o in model TD as free parameters and fitted them to 
the observed quantities. The procedure was iterated by shifting the central colours 
and the colour gradients according to the corresponding uncertainties. We found 
that the change of To and E(B — V)o are correlated, in the sense that (5(logTo) ~ 
— l/3(5(logE(B — V)o), and that the variation of E(B — V) and T from the center to the 
periphery, (5(E(B — V)) and (5(T), are constrained in the ranges 8%<(5(E(B — V))<17% 
and 45%<i5(T)<48% of the central values®. By increasing the inner and outer metallicities 
to Zj = Zo = 3 / 2 Z 0 , the best fit values of E(B — V)i is ~ 0.09, and the following con¬ 
straints are obtained: 30%<i5(E(B — V))<90% and 30%<(5(T)<45%. By increasing r to 
10 Gyr (with Z; = Zo = O.ddZ©), we obtain E(B - V); = 0.26, 51%<(5(E(B - V))<53% 
and 40%<d(T)<45%. To summarize, in order to explain central colours and colour gra¬ 
dients of disks, both age gradients and dust absorption are needed. For a wide range 
of stellar population parameters, the radial age variation is strongly constrained in the 
range [30, 50]% of the central value, while E(B — V)i is greater than 0.1. On the other 
hand, the radial gradient of dust absorption depends critically on the inner metallicity 
and on the time scale of star formation rate. 

So far, we considered only models for which the dust gradient does not change from 
one cluster to the other. By relaxing this assumption, the presence of age gradients is no 
more needed in order to explain all the observed quantities. This is shown by model D: 
a difference of the dust gradient between z < 0.31 and z = 0.64 also fits the data, giving 
central colours (UV - OPT)^ 2 ^ - 2.4, (UV - OPT)^ - 2.5, (OPT - NIR)g - 3.9, 
(OPT —NIR)Qg 4 = 3.8, which are in good agreement (at < 1.5 (t) with the values in 
Table 0 As shown by Peletier et al. 119951 useful tools to investigate the role of dust 
gradients in galaxies are the diagrams showing colour gradients against central colours 
and/or inclination. If dust is the primary effect driving the observed gradients (model 
D), we would expect that redder and more inclined disks have also stronger OPT-NIR 
colour gradients for both AC 118 and EIS0048. The diagrams are shown in Fig. ITTI 
By analyzing the figure, it is important to take into account that the uncertainties on 
central colours and colour gradients are correlated. At least for AC 118, this correlation 
fully explains the tight relation between grad(OPT — NIR) and the OPT — NIR colours, 
while for EIS0048 some correlation could exist. On the other hand, the diagrams with 
inclination vs. grad(OPT — NIR) do not show any correlation at all (cfr. Fig. 1 of Peletier 
et al. ll995|l . These facts seem to suggest that dust is not the primary cause of the observed 
gradients, and that models TD are a more suitable explanations of the present data. 

Finally, some considerations are needed about the ranges of values of 
grad(UV — OPT) and grad(OPT — NIR). We see that age gradients stronger than T1 
can successfully describe the tails toward negative values of these distributions, with the 


® The extrema mark a 10% percentile interval 
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Fig.11. OPT-NIR central colours and logarithm of axis ratio vs. OPT-NIR colour gra¬ 
dients for AC118 (left panels) and EIS0048 (right panels). Dashed lines mark the mean 
values of grad(OPT — NIR). The ellipses in the upper panels correspond to Icr confi¬ 
dence contours and show the correlation of the errors on OPT-NIR colours and colour 
gradients. 

constrain To > 0.5 • Tj set by grad(UV — OPT) at z = 0.64. For what concerns the pos¬ 
itive values of grad(UV — OPT), which correspond to negative OPT — NIR gradients, 
they seem not to be explainable with the models considered above. 

7. Summary and conclusions. 

We have studied the waveband dependence of the shape of the light distribution in cluster 
galaxies from z ~ 0.2 to z ~ 0.64, by using a large wavelength baseline covering UV to 
NIR restframe. New structural parameters have been derived for the galaxy clusters A 209 
at z = 0.21, in UV and OPT restframe (B and R bands), and for EIS0048, at z = 0.64, 
in UV, OPT and NIR restframe (V, I and K bands). Previous determinations of OPT 
and NIR (R and K bands) structural parameters for the cluster AC 118 at z = 0.31 have 
also been used. The present data constitute the first large (N ~ 270) sample of cluster 
galaxies with surface photometry from UV to NIR at intermediate redshifts. The analysis 
has been performed in terms of both structural parameters and colour gradients/central 
colours of galaxies, for the populations of disks and spheroids. 
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7.1. Spheroids. 

On average, the light profile of galaxies is more concentrated at longer wavelengths: 
effective radii decrease from UV to NIR, while Sersic indices increase. The values of 
A(logre) and A(logn) seem not to evolve significantly with redshift, although some dif¬ 
ferences (at ^ 2cr) among the cluster could exist. By averaging the results for A 209, 
AC 118 and EIS0048, we obtain ry^/rOP"^ = 1.2 ± 0.05, n^P'^/n^v = i.Q ± 0.1, and 
r°P'^/rN» = 1.26 ± 0.06, n°P^/nN» = 0.88 ± 0.03. The value of r^P'^/r^iR is fully con¬ 
sistent with that obtained for nearby early-type galaxies by Pahre et al. 119981 who found 
r^/r^ ~ 1.2. The ratio of UV-OPT effective radii at z ^ 0 can be estimated from the 
U- and R-band data of galaxies in rich nearby clusters by Jprgensen et al. and 

amounts to r^/r^ ^ 1.17, which is in full agreement with our results. The bulk of 
spheroids shows negative colour gradients, both in UV - OPT and in OPT - NIR, and 
is therefore characterized by stellar populations redder in the center than in the pe¬ 
riphery. The mean values of colour gradients are grad(UV — OPT) = —0.18 ± 0.04 and 
grad(OPT — NIR) = —0.4 ± 0.1, and are consistent with the values obtained by PDI90, 
grad(U — R) = —0.21 ± 0.02, by Idiart et al. II2()()3II . grad(U — V) = —0.23 ± 0.07, and 
by PVJ90, grad(V — K) = —0.16 ± 0.06, for nearby early-types. Moreover the UV-OPT 
colour gradient is fully consistent with that estimated by Tamura & Ohta Il 2 ()()()ll for 
cluster early-types at z ^ 0.4, grad(B 45 o — I 814 ) = —0.23 ± 0.05. We conclude, therefore, 
that the colour radial variation in spheroids does not change significantly at least out to 
z = 0.64. The comparison of colour gradients and central colours with gradient models of 
stellar population parameters shows that metallicity is the primary cause of the colour 
variation in galaxies, in agreement with previous optical studies at intermediate redshifts. 
The data are well fitted by models with super-solar inner metallicity (Z, = 3/2Zq), and 
with a pure metallicity gradient of about — 0 . 2 dex per decade of radius, in agreement 
with what found by Saglia et al. Il 2 ()()()ll and Idiart et al. onmti . Taking into account 
age-metallicity degeneracies, however, mild age gradients are also consistent with our 
results, the age of stellar populations (referred to z = 0 ) being constrained to vary by 
less than 25% from the center to the outskirts. 

7.2. Disks. 

The light profiles of galaxies become more peaked in the center from UV to NIR. For 
what concerns structural parameters, we find that effective radii decrease at longer wave¬ 
lengths, while Sersic indices become greater. The ratios of UV-OPT parameters and 
OPT-NIR effective radii do not show significant evolution with redshift. The mean val¬ 
ues of rUV/rOPT^ j^uv/jjOPT j.opt/j.nir amount to 1.1 ± 0.035, 0.84 ± 0.03 and 
1.25 ± 0.04, respectively. The most sharp result, however, is the decrease in the ratio 
of OPT to NIR Sersic indices from z = 0.31 to z = 0.64: the value of n®P"^/n^^P varies 
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from ~ 0.7 at z = 0.31 down to ~ 0.45 at z = 0.64, implying that the OPT-NIR structure 
of disks becomes significantly more concentrated at higher redshift. When analyzed in 
terms of colour gradients, the previous results entail a presence of UV-OPT and OPT- 
NIR gradients much stronger with respect to those of spheroids. This is similar to what 
found in studies of colour gradients for nearby galaxies (e.g. Peletier et al. rrn^mii . weUnd 
that the UV-OPT gradient does not vary significantly with redshift, and amounts to 
^ —0.42, while, on the contrary, the value of grad(OPT — NIR) strongly decreases from 
~ —0.7 at z = 0.31 down to ~ —1.2 at z = 0.64. Moreover, disks have redder colours in 
the center, comparable to those of spheroids. We analyze these results by using simple 
gradient models describing the effect of variations of stellar population parameters as 
well as dust extinction throughout the galaxies. It turns out that two different kinds of 
models fit the present data. 

1. Age gradient models with (i) a protracted star formation rate both in the inner and 

in the outer regions, and with (ii) a significant amount of dust extinction in the 
center. The property (i) keeps the variation of grad(UV — OPT) mild and allows the 
observed evolution of grad(OPT — NIR) to be reproduced, while the property (ii) is 
needed in order to describe the red central colours in disks. The better constrained 
quantities are the age gradient, in the range [30,50]% of the central value, and the 
inner extinction, E(B — V)i>0.1. On the other hand, the amount of dust in the outer 
region depends critically on the value of other stellar population parameters (inner 
metallicity and time scale of star formation). It is interesting to note that the age 
gradient model predicts an UV — OPT colour gradients of ^ —0.35 at z ~ 0, which 
is in very good agreement with that measured by Gadotti & dos Anjos llWll for 
a sample of nearby late-type spirals, grad(UV — OPT) ~ —0.33 (see their Table 2). 
The existence of age gradients in disk dominated galaxies is in agreement with what 
found by de Jong for nearby spirals, while the fact that dust absorption plays 

an important role in disks is in agreement with the findings of Peletier et al. II 0051 
for Sb and Sc galaxies at z ~ 0. Moreover, age gradients are a robust prediction of 
hierarchical models, in which disks are accreted by the condensation of gas from the 
surrounding halos. 

2. Pure extinction models, in which the gradient of E(B — V) increases by a factor ~ 2 
for EIS 0048 with respect to the other clusters. 

In order to discriminate between the two models, we analyze the relations between colour 
gradients and central colours/axis ratios. The present data indicate that galaxies with 
higher inclination do not have stronger colour gradients, as one would expect if dust 
gradients would be the main effect driving the observed gradients. Although some cor¬ 
relation between colour gradients and inclination could be hidden by the uncertainties 
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Fig. A.l. Image of the central region of EIS 0048 in the high resolution (left panel) and 
in the standard resolution (right panel) I-band images. The images are normalized to the 
same intensity scale. Note the spatial scale in the lower-left panel. 

on the above quantities, we argue that age gradient models could give a more suitable 
explanation of the present data. 

Since hierarchical merging scenarios predict strong differences in the properties of 
field and cluster galaxies, it will be very interesting to analyze the role that environment 
has on the internal colour distribution of galaxies. Although we show by a statistical 
subtraction procedure that our results are not affected from field contamination, we 
cannot discriminate between the properties of field and cluster galaxies. In order to 
address this point, spectroscopic data will be needed. 
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Appendix A: Reduction of FORS2 I-band images. 

New I-band images of the cluster of galaxies EIS 0048 have been collected at the ESO Very 
Large Telescope (VLT) during August 2001 with the FORS2 instrument, by using the 
high resolution observing mode (pixel scale 0.1'Ypxl). The nights were photometric with 
excellent seeing conditions (0.3" < FWHM < 0.4"). The data consist of three pointings, 
for each of which we obtained hve dithered exposures of 360 sec each. Data reduction 
was performed as described in LMI03, by using IRAF and Fortran routines developed 
by the authors. Images were bias subtracted and corrected for flat-field by using twilight 
sky exposures. The frames of each night were median combined in order to obtain a 
super-flat frame, that was used to improve the accuracy of the flat correction (at better 
than ~ 0.5%) and to fully remove the fringing pattern in the images (see LMI03 for 
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details). The frames were registered by integer shifts and combined by using the IRAF 
task IMCOMBINE. The photometric calibration was performed into the Johnson-Kron- 
Cousins system by using comparison standard fields from Landolt lITM^ observed during 
each night. The accuracy on the photometric calibration amounts to ~ 0.005 mag. In 
Fig.lO we show, as example of the final quality of the image, the central region of the 
cluster. The image is compared with that obtained in the I band under more ordinary 
seeing conditions (FWHM ~ 0.8") with the standard resolution observing mode (pixel 
scale of 0.2"/pxl, see LMI03). 

Appendix B: PSF of the FORS2 l-band images. 

The PSF of the FORS2 images was modelled by using a multi-Gaussian expansion, taking 
into account the effect of pixel convolution (see LBM02). Stars were selected according 
to the stellar index (SG) of SExtractor (Bertin & Arnout flOOfill . We considered only 
objects with / < 20 and SG > 0.95, in order to obtain a reasonable accuracy in the 
fitting and to minimize the contamination by extended sources. In Fig. IH3I we show 
the contour plot for a star in the cluster field and the corresponding circular model, 
obtained as described in LBM02. It is evident that the star is asymmetric in a direction 
at ~ —45 degrees with respect to the horizontal axis. This effect can be explained by the 
fact that five dithered exposures were available for each pointing and the seeing was not 
exactly the same for each of them. Since the images were combined by integer shifts, the 
seeing difference can easily account for the PSF asymmetry. We found that this effect can 
introduce a significant bias (10- 15%) in the measurement of the structural parameters 
and, therefore, we considered more suitable models of the PSF, by using the rotation 
angle of each Gaussian as a free parameter in the fit and by introducing another free 
parameter to describe the deviations of the PSF isophotes from the circular shape. In 
Fig. EH we see that the asymmetric model matches very well the star isophotes. For 
each image, we fitted first each star separately, in order to investigate possible variations 
of the PSF across the chip, and then, since no significant dependence on the position was 
found, we constructed the PSF model by fitting all the stars simultaneously. 

Appendix C: Surface brightness profiles. 

The UV, OPT and NIR surface brightness profiles of the N = 273 galaxies of A 209, 
AG 118 and EIS0048 are shown in Figs. G.l, G.2, and G.3^°, respectively. The ID 
profiles are plotted down to a signal to noise ratio S/N = 1. For each galaxy, we sampled 
the observed, fitted, and de-convolved images along ellipses with center coordinates, axis 
ratio, and position angle given by the two-dimensional fit, and derived the mean surface 
brightness ^(r) as a function of the equivalent radius r of each ellipse. We considered only 


Figures are available in electronic form at http://www.edpsciences.org 
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Fig. B.l. Modelling of the I-band PSF for EIS0048. The contour plots refer to a star 
in the cluster field (continuous lines) and to the PSF models (dashed curves). Left and 
right panels show the symmetric and asymmetric models, respectively. 

the pixels not masked in the 2D fit and excluded ellipses with a low fraction (< 30%) 
of unmasked pixels. For the observed profiles, the error bars on /i(r) were computed 
by taking into account the surface brightness fluctuations along the ellipses as well as 
the local background noise of each image. The profiles of galaxies belonging to multiple 
systems were obtained by first subtracting the 2D models of the companions from the 
observed galaxy image. 
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